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A density matrix theory of electron transport and optical gain in quantum cascade 
lasers in an external magnetic field is formulated. Starting from the general quantum 
kinetic treatment, we describe the intra- and inter-period electron dynamics at the 
non-Mar kovian, Markovian and Boltzmann approximation levels. Interactions of 
electrons with longitudinal optical phonons and classical light field are included in the 
present description. The non-Markovian calculation for a prototype structure reveals 
significantly different gain spectra in terms of linewidth and additional polaronic 
features in comparison to the Markovian and Boltzmann ones. Despite strongly 
opposed interpretations of the origin of the transport processes in the non-Markovian 
or Markovian and the Boltzmann approaches, they yield comparable values of the 
current densities. 

I. INTRODUCTION 

The rapid experimental progress in the field of quantum cascade lasersi (QCL's) has ini- 
tiated considerable theoretical activity to explain the underlying physical phenomena and 
improve their performance by design optimization. To date, both semiclassical and quantum- 
mechanical theories of carrier transport in QCL's without magnetic field have been proposed. 
Semiclassical ones are based on the assumption that coherent processes in QCLs are negligi- 
ble and electron transport occurs via scattering processes only. They rely on the Boltzmann 
transport equation, for solution of which a few approaches may be employed. The Monte 
Carlo method is a stochastic approach which simulates the trajectories of a representative en- 
semble of carriers .-i^i^ An assumption that the carrier distribution in any particular subband 
can be approximated by Fermi-Dirac statistics allows the Boltzmann equations to be replaced 
by simpler and less computationally demanding rate equations.-i^ The quantum-mechanical 
models enable the description of phase coherence as well as incoherent scattering processes, 
and they have been formulated using the density-matrix or nonequilibrium Green's function 
approach.-iSi'^ii^iiiii^ii^ii^ Comparison of the results obtained with the Boltzmann and density 
matrix approaches in a mid-infrared QCL performed by lotti at al.- showed that quantum 
corrections to the current density are negligible. However, the analysis of gain spectra in the 
nonequilibrium Green's function description demonstrated that the negligence of coherences 
between QCL states results in significantly broader linewidths-^^ Also, it has been argued 
that coherences are not irrelevant for transport in terahertz (THz) QCL's where the small 
anticrossing energies allow for resonant tunneling.— Moreover, a very recent study^^ gave 
an interpretation where the current across QCL's is entirely coherent. In addition to the 
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studies concentrated on QCL's, there is a mounting theoretical evidence for the presence 
of quantum coherence features in hnear absorption spectra and nonhnear ultrafast optical 
response for intersubband transitions in unbiased quantum wells (QW's).— '1'^'^^'^^ 

Furthermore, experimental interest in the QCL performance in a magnetic field stimu- 
lated theoretical efforts to describe the influence of magnetic field on the physical processes 
involved. However, since this research topic emerged recently, very few theoretical studies 
of QCL's in a magnetic field, compared to the amount of those for QCL's without magnetic 
field, have been reported. Most of them were focused on the modeling of various scatter- 
ing rates (electron-longitudinal optical phonon,— i2ii2^i2^ electron-electron,— and interface 
roughness^) and the calculation of these scattering rates between the upper and the lower 
laser levels. Modeling of the active region of QCL's, including electron-longitudinal optical 
(LO) phonon and electron-longitudinal acoustic (LA) phonon scattering, and assuming a 
unity injection approximation, has also been reported.— Finally, a semiclassical model 
of the electron transport in a magnetic field based on the Boltzmann equation has been 
proposed.— Apart from the work done on QCL's in a magnetic field, a few theoretical in- 
vestigations of QW systems subjected to a magnetic field, based on both the density matrix 
and nonequilibrium Green's function approaches, have been reported, and confirmed the 
importance of quantum coherence effects on the ultrafast time scales.— 

Currently, no experimental or theoretical data on coherent phenomena in QCL's in a 
magnetic field are available. Since the energy spectra in such structures is discrete, it is rea- 
sonable to expect that coherent effects are more significant than for QCL's without magnetic 
field. The aim of this work is to present a quantum-mechanical theory of transport and gain 
properties of QCL's in an external magnetic field. For that purpose, we derived quantum 
kinetics equations for QC structures in a magnetic field, based on the density matrix formal- 
ism, which include interaction of electrons with LO phonons and optical field. Furthermore, 
we obtained the corresponding equations in the Markovian approximation, from which the 
semiclassical Boltzmann transport equations can be recovered. A comprehensive analysis 
is performed for an example GaAs/Alo.aGao.yAs QCL and nonequilibrium steady state re- 
sults obtained from all three approaches (quantum kinetic, Markovian and Boltzmann) are 
compared. 



We consider electrons in the conduction band of a QCL in a magnetic field applied 
in the direction perpendicular to QW layers {z axis). Such magnetic field splits the in- 
plane continuum of quantized subbands into Landau levels (LL's), additionally described by 
Landau and spin indices.— Within the effective mass and envelope function approximations, 
and neglecting the spin splitting, the energy of the jjth LL originated from the mjth state 
(subband), in further considerations denoted with a shorthand subscript i, i = \mi,ji), reads 



where is the energy of state mj, h is the reduced Planck's constant, e is the electron 
charge, B is the applied magnetic field, and m* is the electron effective mass (taken to be 



II. THEORETICAL CONSIDERATIONS 



A. Quantum kinetic equations 




(1) 
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equal to 0.067 in free electron mass units). For the magnetic vector potential A given in the 
Landau gauge (A = Bxey), the envelope wave function of the ith LL takes the form 

^,,fe(r) = u,^ {x - f{k))^^Xz)^, (2) 

where k is the wave vector of the electron, uj. {x — f{k)) is the wave function of the harmonic 
oscillator with f{k) = k/ (eB/h), iprmiz) is the wave function of the m^th size-quantized state, 
and Ly is the dimension of the structure along the y axis. 

The model Hamiltonian of the system described above reads 

H = Ho + Hel + Hep. (3) 

The first term represents the Hamiltonian of noninteracting electrons and phonons in ap- 
plied electric and magnetic fields. The second and the third term describe electron-light 
and electron-LO phonon interactions, respectively. In this first step towards formulating a 
density matrix theory of QCL's in a magnetic field, we do not consider other interaction 
mechanisms of electrons (with LA phonons, ionized impurities, interface defects, and other 
electrons). The Hamiltonian of free electrons and phonons reads 

Ho = J2Eicy.,k + J2Ej>% (4) 

i,k q 

where Ei is the energy of the ith electron state, i?q is the energy of the phonon of a wave 
vector q, and clj^ (6jj) and Ci^k (^q) represent creation and annihilation operators of the 
electron (phonon), respectively. Electron-light interaction in the dipole approximation is 
given by the following Hamiltonian 

Hel = E eA^Vlfcy,,,,, (5) 

i,j,k,k' 

where A/j represents the magnetic vector potential of a monochromatic light wave incident 
on a QW structure, given in the Coulomb gauge, and the velocity matrix element is found 
according to 

\1f= Jdr^l,{r)vo^,,,ir). (6) 
The velocity operator for the nonilluminated system vq may be represented as 

1 ^ eA 

vo = — P + — , 7 

m m 

where p is the momentum operator. The Hamiltonian describing electron-phonon interaction 
can be cast in the form~'^ 

Hep= E +£'fc!q,fe'cifc'^i^*.fc)' (S) 

«j,fc,fc',q 

with 

9t^,k' = 9c,j dv%^,{r)e^^-^^,,,.{v), (9) 
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where represents the couphng factor. For electron-LO phonon interaction, the couphng 
factor reads 

hujj 



9ci = -ie 



LO / -1 -1 



2V 



Coo - 



1/2 ^ 



(10) 



where the energy of each phonon mode is considered to be approximately constant {hui^o), 
V is the volume, and and are high-frequency and static permittivity, respectively. In 
the case of QW's in a magnetic field, the phonon coupling factor gV ^, may be written as 



9Lk' = 9jdv ul {x - f{k)) '-^r^^{z)e^^'^-^^'^yy^'^^^\,^ {x - f{k')) ^^„,.(^) 



-i{k-qy-k')y 

9c, I dy / dx ul (x - f{k)) e'^-^u,^ (x - f{k')) I dz rn^MY'-'^m.^z) 



Ly 



9n^k',k-qyHj^j^{k, k', q^)G 

(11) 



where Hj-j.{k, k', qx) = J dx ul (x — f{k)) e^^'^^Uj^ (x — f{k')) is the lateral overlap integral 
and G„nmj{qz) = J dz ip^_{z)e'^'^''^ilJmj{z) is the form factor. 

In the density matrix approach, single particle density matrices like the intraband elec- 
tron density matrices fi-^i2,k = (c|i,A:Cj2,fc) or the phonon occupation number nq = (blfiq^ 
represent fundamental physical quantities. Their diagonal elements determine the occupa- 
tion probabilities of the states, while the nondiagonal elements correspond to the electron 
polarizations between two states, and are related to the property of quantum-mechanical co- 
herence (superposition). In this work, a thermal equilibrium of phonons is assumed, hence 
equations of motion for electron density matrices are sufficient for the description of the 
system. 

In the derivation of the time evolution of single particle density matrices, one starts 
with the Heisenberg equation of motion.— i^i The time evolution due to the Hamiltonian of 
noninteracting electrons and phonons is given as^"^ 



d 



dt' 
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i,fili2,k Ho ~ -fc (-^«2 -^ii)/iii2,fc (12) 



ih 



The equation of motion in the case of interaction of electrons with light polarized in the 
z-direction reads^^i^i^ 



d e 

fhi2,k H^j^ — TT ^i? (^2«3/n«3,fc ~ ^3*1 /«3«2,fc) 5 (l^) 

«3 



dt 

with the ^-component of the velocity matrix element 



I 

Viii2 — ^("^^i ~ -^i2)diii2^ji^,ji^^ (14) 

where di^i^ is the matrix element of the operator of the ^-coordinate. In the present analysis 
we do not consider any optical-cavity effect and look for nonequilibrium steady state popu- 
lations and polarizations when = 0. However, electron-light interaction is essential for 
the calculation of optical gain. 

In the quantum kinetics equations for electron-phonon interaction, phonon-assisted ma- 
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trices, given with expectation values of 3 operators s^^q^,, = \c\^ JjqCi2,k'jj appear, which 
correlate an initial state consisting of one electron in the state i2ik' and a phonon with a 
wave vector q to a final state with only one electron in the state ii,kM^ Furthermore, 
the temporal evolution for the phonon-assisted matrices involves expectation values of four 
operators, and so on. The resulting infinite hierarchy of equations needs to be truncated in 
order to access the problem numerically. The first order contribution, obtained by neglect- 
ing all correlations between electrons and phonons in the spirit of the correlation expansion 
approacb^^-^ (4^^%, ^ (4,fcCi2,fc) (V.) ^fc'.fc^g^.o = /ni2,fc5q..^fc',fc5g„o), vanishes if a ther- 
mal equilibrium of phonons is assumed.— The next order in the hierarchy is obtained by 
taking into account deviations of the phonon-assisted density matrices from the first order 
factorization ^s'^.'^q^,, = s^^q^,, — fi^i2,kBci.^.^5k',k5qyfl- Then, the following equations for ^s^'^q j^., 
are obtainecP^'^ 

+ ^ X! dk^aj:' [(""-O + ^) fills, k{^U,i2 ~ fiii2,k') " '^o/j4i2,fc' ('^ii.is ~ fhi5,k)] , 



d_ 



f . \ " /n^a's xJih _|_ ^«3«2* Xoian* _ ^«3H Ae*3«2 _ „n«3* x„«2i3* ^ 

Jhi2,k Hep ~ ^f. \^yfc,q,fc'"'^fc,q,fc' i/fc'.q.fc^-^ifc'.q.fc y k' ,(i,k'^ ^ k' ,q,k i/fc,q,fe'"'^A:,q,fc7 ' K^"^) 



ih . 

«3,K ,q 



where hq denotes the equilibrium phonon density given by the Bose-Einstein factor. The 
terms in the equation for 5s^^q ^/ are due to the Hamiltonian of free electrons, and of electron- 
LO phonon interaction, respectively. The equations for phonon-assisted matrices should, in 
principle, contain a term which describes their time evolution due to electron-light interac- 
tion, here relevant only for the calculation of linear optical gain. However, the coupling of 
the light field to the phonon-assisted matrices in QW's is a higher-order effect^S and may 
be neglected.— Its inclusion for complex structures like QCL's in a magnetic field would 
result in a computationally inaccessible task.— 

Insertion of higher order terms in the equations for phonon-assisted density matrices 
should be performed in a self-consistent manner,— however, in the system considered, with 
several subbands and LL's originating from them in each period of the cascade, this would 
be extremely computationally involved.—"^ Conversely, discarding these effects leads to 
numerical instabilities in the actual computation. Therefore, a phenomenological damping 
constant 'jph was introduced, representing higher order correlations.— We have verified that 
the convergence of our results may be achieved for sufficiently large values of 'jph (~ 1 meV). 

It is shown in Appendix A that, in the present description, /j^jj^fc is constant for all values 
of the wave vector k, and may be expressed as /iiij.fc = tt_B^nj25 where = nh/eB, and 
^4142 = J2k' fhi2,k' /LxLy. The diagonal element ria represents the electron sheet density in 
the ith LL. From the derivation given in Appendix A, quantum kinetics equations including 
all the aforementioned interactions amount to 

d 1 e 

~n^iii2 '^i.Bi2 Ei^Tli-^i^ H — 2.-^R (^2i3'^n«3 ^1311^^12) 

dt in in 

«3 

~^jfl S {^^ihUib^ -^hhUis ~^ ^13121514^-^13^11514 " ^«3«i«5«4'^-^«3«2«5M ~ ^11131415^-^12131415} ' 

i3,M,«5 
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with the quantities ^i^j^jj^is associated with the phonon-assisted matrices through 

^^k,q,k-qy = X! 9qHj,_^{k, k — Qy, <ix)G rrii rni^{<iz)5Ki^i^i^i^. (17) 



14, Jg 



The quantum-kinetic dynamics is essentially non-Markovian, since the time evolution of 
density matrix elements depends on their values at earlier times i.e. on the memory of the 
system. 

From discussion in Appendix A it follows that the time evolution of electron populations 
and polarizations in the Markovian approximation may take the form 

d 1 , „ „ s e 



-n 



dt in in 

«3 



T^OUt 

■pin _ 



Terms r°"*''/J^g have similar form as scattering rates in the Boltzmann approach, and hence 
may be referred to as generalized out /in scattering rates. The Markovian approximation 
neglects the memory time of a scattering process, which is related to energy-time uncer- 
tainty.—'^'^ Scattering and dephasing processes are then restricted only to energy conserving 
transitions between single-particle states. For discrete energy spectra in QW's in a magnetic 
field, the electron-LO phonon interaction is thus almost fully suppressed, if broadening is not 
taken into account. Therefore, a Lorentzian with the full width at half maximum (FWHM) 
of 7 ~ 1 meV was used to model the LL broadening in the Markovian description. In the 
semiclassical limit, which may be obtained by neglecting nondiagonal matrix elements,—*^ 



131415 

Si-Ei^ + huj-^o + EijWi^i^i^i.^{no + 1) + 5(-^i5 - hui^o + EijW*^i^i^i^nQ\ , 
6{-Ei^ + hui^o + E,,)Wi,,,^,i,no + 5{-Ei, - hui^o + ^44)^^*3^1*544(^0 + 1)1(.18) 



the derived Markovian equations reduce to the Boltzmann equations given in Ref. |29 



Due to the periodicity of the QCL structure, its energy states are invariant upon trans- 
lation per potential drop across a period, while the wave functions are invariant upon 
translation per period length. Therefore, each period has an identical set of LL's, 
with identical density matrix elements {rii^i^ = n(^i-^+kN){i2+kN), k = 0, ±1, ±2, ...), and 
phonon-assisted matrices {5Ki^i^i^i^ = SK(^i^+kN)ii2+kN)ii3+kN)iu+kN))- This also accounts 
for the quantities characterizing scattering processes {Ti^i^i^i^ = T^i^+kN)(i2+kN)(i3+kN){i4+kN), 
^ni2i3M = ^{ii+kN)(i2+kN){i3+kN){i4+kN)) and the velocity operator {Vi^i^ = V(^i^j^kN){i2+kN)) ■ 
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Since the wave functions are well localized within their periods, the tight-binding descrip- 
tion may be introduced, by accounting for the interaction between the nearest neighboring 
periods only. Hence, we consider the density matrix elements which couple LL's within one 
period, as well as the elements which couple those LL's with LL's belonging to the nearest 
neighboring periods. Also, we take into account those quantities Wij^ ^isiiis Fj^ jjgj^jg with 
the properties zi,2 — *3 = 0, ±1 and — = 0, ±1, see Eq. ( IT6il . and write Eqs. ( fT6l) and 
( ITSll for all possible combinations of indices ii — which satisfy these conditions. After 
exploiting the property of shift-invariance of all the aforementioned quantities, the system 
of equations in the quantum- kinetic/Markovian approach may be reduced to contain only 
the density matrix elements of interest. Again, the Boltzmann expressions may be recovered 
from the Markovian ones. 

In the quantum-kinetic and Markovian representations, the number of the density matrix 
elements to be calculated is of the order of A^^, and the number of the quantities associated 
with phonon assisted matrices (the quantum kinetics case only) and scattering rates is of 
the order of A^"^. Obviously, the calculation of population and polarization dynamics for 
the QCL's with many energy states and LL's stemming from them is extremely challenging. 
Therefore, in our analysis, we restrict to the case of a QCL with a small number of energy 
levels per period, and subjected to relatively large magnetic fields, characterized by a small 
number of LL's stemming from those levels which are relevant for the transport. Here 
we took 10 lowest LL indices, after checking that this number of LL's is sufficient for the 
considered structure. 

Stationary solution of Eqs. (fTB]l and ffTS]) is found by tracking their time evolution, starting 
from an initial condition that all electrons are in the fundamental ground state LL (and hence 
all the polarizations and phonon-assisted matrices are equal to zero), and integrating in time 
until the steady state is reached. This method proved to be extremely reliable in terms of 
convergence for solving large systems of nonlinear equations, in contrast to gradient-based 
methods. The integration is performed by using a Runge-Kutta method with adaptive step 
size control, which considerably speeds up the process. 

Since the quantities associated to scattering processes Wi^i^i^i^ and rj^j^jgi^ are different 
from zero only if the condition jj^ = ji2 + Jm is fulfilled, our choice of an initial condition 
leads to the steady-state solution in which the polarizations nj^jj are not equal to zero only if 
jii = ji2- Although this result may, at first, seem to be a peculiarity of the initial condition, 
it can also be regarded as a solution of the reduced description of the systems of Eqs. ffTB]) 
or (fT8l) . which includes only those density matrix elements rii^i^ with ji^ = ji^. Careful 
examination of Eqs. (fT6|) and (fTS!) suggests that, if the terms such as asni^is, iu 7^ Jia; 
are much smaller than 1 (i.e. 5i^i.^ for = i^), the terms such as n^^ig, ji^ ^ ji^, do not 
infiuence the time evolution of the elements rii-^^i^, ji^ = ji^- Therefore, in the lowest order 
approximation, the analysis of the density matrix dynamics may be limited to the terms 

"^1112 ) Jil Ji2 ■ 



The current density may be estimated from the expectation value of the carrier drift 
velocity 



B. Current density 




(19) 
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In the density matrix formalism, the drift velocity may be calculated according to^ 



i\,i2,k 



(20) 



where Vi^i^ is the drift velocity matrix element given with 



The drift velocity density matrix elements then reads 



I 

r 



12 I ^^-^R^ii,i2- 

m* 



(21) 



(22) 



The first term is due to the Hamiltonian of noninteracting electrons, while the second one 
is due to the electron- light interaction. Since the scattering potential of any interaction 
(electron-phonon, electron-electron etc.) depends only on the position r, but not on the 
momentum p, it follows that [Hep, z\ = 0, and its contribution to the drift velocity vanishes. 
Starting from Eqs. (fT9!) and (!20!) . and using the assumption that only the density matrix 
elements between LL's within a period or LL's localized in that period and its nearest 
neighbors are not zero, the current density finally may be written in the form 



, e ^ r 



(23) 



n,«2=l 



where d is the length of a period. The Boltzmann expression for the current density may 
be derived from Eq. fl251) by representing the non-diagonal density matrix elements, in the 
first order approximation, in terms of the diagonal ones^i 



J 



d 



N 



{zf - Zi) [niWif{l - asUf) - nfWfi{l - a^ni)] 



,/=i 
Ui</) 



TV 



+ {zf+N - Zi) \niWi(f+N){l - aBUf) - nfW(f+N)i{l - asrii 



(24) 



where rij = nu represents the population of the ith LL and Zi = da = {i\ z 



C. Gain spectra 



The gain spectra in the quantum-kinetic (non-Markovian) description and in the Marko- 
vian approximation may be estimated from the linear response of nonequilibrium stationary 
populations and polarizations to a small optical perturbation. The relationship between the 
linear variations in the polarization due to the applied optical field AP(a;) and the current 
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density AJ{uj) gives the following expression for the susceptibilityii 

_ APjoo) _ i AJjio) 

The gain coefficient then may be found fron>2^ 

iulm[xiuj)] 



(25) 



c n 



(26) 



where n is the refractive index of the system material. 

If Fourier transform of the electric field of light is given with 



E(t) = e. 



duj 
2^ 



E{ij)e 



-iuit 



(27) 



then Fourier transform of the corresponding magnetic vector potential in the Coulomb gauge 
is represented as 



A«(t)=e.y- 



duj E{ijj) 



-iu)t 



lUJ 



(28) 



The linear changes of intra-period elements in the frequency domain Ani^i^{uo), for the 
quantum kinetics case, may be written as 



- iujAni.iJuj) 



is 



-UJ 



*3,*4,«5 



, . 1 



Ei^ + huji^o - Ei^ -hu - ih-fph 
{{no + 1) [An^,i^{uj)6i,^,^ - aB{nli^Ani.^i^{uj) + nl^^Am.i^iu)) ^ 

-no [Ani^i^{uj)6i,^i^ - aB{nl-^Ani^i^{uj) + n° An,^i,(a;))] } , 



(29) 



where n° represents the steady-state value of the density matrix element between LL's 
ii and 22- Similar equations which include all possible combinations of ii — % as discussed 
in Subsection III Al also need to be taken into account. In the Markovian case, taking into 
consideration the property of the density matrix elements that ni^i^{t) = n*_^i^{t), the equa- 
tions of motion of intra-period density matrix elements transform to the frequency domain 
according to 



iuArii^i^^u) 



is 



«2«3"'liJ3 '^«3n"j3J2y 



+ E [-r°2t^4i5(^^n*5(^)'^^4,i3 - aBK^5^'^^4^3(^) + ^^U^S (^) ) ) 



«3,*4,i5 



■^out * 



{Ani^i^{Lj)6i,^i^ - asin^i Aui^i^iu) + nl^ Am^i^iu))) 
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(30) 



Current density in the frequency domain in both the non-Markovian and Markovian descrip- 
tions may be obtained from 



e \JL 1 



N 



n=i 



«1,«2=1 



In the Boltzmann description, the gain coefficient may be given as 

2 N 
TIC r 

9{^) = E (E, - Effsgn{E, - Ef)dl^^^^5{\E, - Ef\ - hu) 

fsgn{E, - Ef+N)d^ 



(31) 



neocioh d ^ 

+ {Ei - Ef+NYsgn{E, - Ef+N)dtn^,mf+J{\Ei - ^/+7v| - ^^) 

+ (E,+;v - E^)\gYy{E,^^ - i?/X^^,^/(|i?.+iv - Ej\ - hu) 



(32) 



The deha function in the gain coefficient expression is modeled by a Lorentzian with the 
same FWHM as for electron-LO scattering rates. 

III. NUMERICAL RESULTS AND DISCUSSION 



As a prototypical system, we consider a QCL design which comprises a three-level scheme, 
and employs LO-phonon depopulation of the lower laser level to the ground state. No injector 
region is present and efficient injection into the upper laser level is enabled by its alignment 
with the ground level of the preceding period. The QCL period consists of two QW's (see 
Fig. [T]), one of which confines the ground and lower laser levels, whose energy difference is 
set to be approximately one LO phonon energy. The upper laser level is localized in the 
other well. This structure was chosen to be examined, instead of existing QCL's already 
investigated in the presence of a magnetic field, due to its simplicity, and the dominant 
influence of electron-LO phonon interaction on the electron population dynamics, as will be 
explained in what follows. 

The conduction band proflle and electronic structure of the QCL in zero magnetic fleld 
and an electric fleld of 16.2 kV/cm, is given in Fig. [T] One QCL period includes a 2.8 nm 
Alo.3Gao.7As barrier, followed by a 9 nm GaAs well, a 1.4 nm Alo.3Gao.7As barrier and a 
17.4 nm GaAs well. States 1, 2 and 3 represent the ground, lower laser and upper laser levels, 
respectively, and states 1' and 3" represent the ground level of the preceding period and the 
upper laser level of the following period. The doping density was chosen to be 10^^ cm~^, in 
order to achieve relatively high gain in the THz range, while still having a small influence on 
the effective conduction band potential and making electron-electron processes less relevant. 
The temperature was set to 4 K, since QCL's in a magnetic fleld are usually operated at low 
temperatures. The transition energy between the upper and lower laser levels is 15.2 meV, 
and the energy difference between the ground state and the upper laser level of the following 
period is 2.6 meV. 
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T ' 1 ' r 




z [nm] 

FIG. 1: A schematic diagram of the conduction band profile, size-quantized energy levels from 
which Landau levels originate and squared wave functions for one full period and parts of adjacent 
periods of the GaAs/AlGaAs QCL for zero magnetic field and an electric field of 16.2 kV/cm. 
States 1 and 1' (solid line), 2 (dash-dotted line), 3 and 3" (dashed line) denote the ground, lower 
laser and upper laser levels, respectively. State 1' belongs to the preceding period, while state 3" 
belongs to the following period. 

Relatively strong electron-electron scattering occurs between the ground state and the 
upper laser level of the next period, due to a large overlap and small energy difference. 
However, the semiclassical calculation performed for a similar structure showed that the 
electron-LO phonon scattering rates from the lower laser level to these levels, also relevant 
for the distribution of electrons between them, are considerably larger.— Electron-electron 
processes between other states are less important, due to the large energy spacing of sig- 
nificantly populated LL's. Consequently, the population dynamics are not significantly 
influenced by electron-electron scattering. We make an assumption that the same accounts 
for the polarization dynamics and neglect electron-electron interaction hereafter. 

A. Electron populations 

The populations of all LL's associated with the ground state and the upper and lower 
laser levels, calculated using the non-Markovian, Markovian and Boltzmann model of elec- 
tron transport, as functions of magnetic field, are shown in Fig. [2l In the non-Markovian 
treatment, we used the values of the damping parameter of Ti-^ph = 1 meV, and fi^^ph = 2 meV, 
while in the Markovian and Boltzmann description, the Lorentzian FWHM of 7 = 3 meV 
was taken. Regardless of the model used, their dependences on the magnetic field are gener- 
ally similar. For certain magnetic fields (4.6 T, 9.2 T), the energy difference between some 
LL's stemming from the upper laser level and the ground state becomes equal to one LO 
phonon energy, thus the electron-LO phonon interaction between them increases consider- 
ably. Consequently, the population of all LL's stemming from the upper laser level decreases 
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reaching its minimum, while the opposite happens to the population of the LL's stemming 
from the ground state. Conversely, for intermediate magnetic fields (6.2 T), the population 
of the upper laser level is increased, while the ground state is depopulated. The population 
of the lower laser level practically does not change with magnetic field. 




FIG. 2: (Color online). The electron population over QCL states (all Landau levels) vs magnetic 
field. States 1, 2, and 3 represent the ground state, the lower laser level, and the upper laser level, 
respectively. Ns is the total sheet density of electrons per period. Solid, dash-dotted, dashed and 
dash-double dotted lines represent non-Markovian {h'jph = 1 meV, and h'jp^ = 2 meV), Markovian 
(7 = 3 meV) and Boltzmann (7 = 3 meV) results, respectively. 

The populations obtained from the Markovian and Boltzmann description do not differ 
much, except in the range of magnetic fields between 6.2 T and 8 T, see Fig. [2l In the 
fully nondiagonal Markovian (and non-Markovian) approach employed here, the coupling 
between populations and polarizations is accounted for, which results in the presence of phase 
coherence in the stationary state. In other words, during the time evolution of the system, 
scattering processes from populations to polarizations create a certain amount of polarization 
in the steady-state conditions. These and reversed processes (dephasing from polarizations 
to populations) may have an apparent impact on electron populations. More pronounced 
differences between the Markovian and Boltzmann populations, for magnetic fields between 
6.2 T and 8 T, indicate a stronger interplay between populations and polarizations in the 
Markovian description and, hence, larger values of polarizations. Indeed, Figs. [Hand [5] show 
that the polarizations between all LL's stemming from any two laser states are increased for 
these magnetic fields. 

On the other hand, the populations calculated from the quantum-kinetic model, for both 
values of the damping rate, show significant departure from the Markovian or Boltzmann 
results. For h'jph = 1 meV {hjph = 2 meV), the amplitude of the oscillations of the upper 
laser level and the ground state populations is considerably larger (smaller) than in the 
Markovian or Boltzmann case, see Fig. [21 Such a drastic contrast between the results of the 
quantum-kinetic simulation for different damping values reflects the fact that these quan- 
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titles actually describe colllslonal broadening,— to which the LL populations are extremely 
sensitive. From that point of view, the difference In the population dynamics between the 
quantum-kinetic treatment and the Markovlan or Boltzmann approach for 7 = 8 meV may 
be explained In terms of correspondence of h'^ph = 1 meV {Ti'-jph = 2 meV) to a smaller 
(larger) broadening, and hence, a smaller (larger) Lorentzlan FWHM. Indeed, the Marko- 
vlan and Boltzmann calculation for different values of FWHM showed that the LL popu- 
lations for 7 = 2 meV and 7 = 4 meV are very similar to the quantum-kinetic results for 
Ti^ph = 1 meV and Ti'-^ph = 2 meV, respectively, see Fig. [31 Nevertheless, this does not mean 
that there Is one-to-one correspondence between these electron transport models for such 
values of FWHM and damping parameter. The non-Markovlan description accounts for the 
memory of scattering and dephaslng processes, which corresponds to a quantum-mechanical 
energy-time uncertainty. In comparison, the Markovlan or Boltzmann dynamics take Into 
account only energy-conserving processes, here, however, relaxed by the assumption that 
all LL's are broadened. This may lead to different values of polarizations, gain and cur- 
rent although the populations are quite similar. In order to give a fair comparison between 
different models presented here, we restrict our further analysis to a choice of the pairs of 
phenomenologlcal parameters for which the populations are almost Identical (7 = 2 meV 
and fi'~^ph = 1 meV; 7 = 4 meV and Ti'^ph = 2 meV), and then we compare the results for 
other physical quantities. 
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FIG. 3: (Color online). The electron population over QCL states (all Landau levels) vs magnetic 
field. States 1, 2, and 3 represent the ground state, the lower laser level, and the upper laser 
level, respectively. Ns is the total sheet density of electrons per period. Left: Solid, dashed and 
dash-double dotted lines represent non-Markovian {hjph = 1 meV), Markovlan (7 = 2 meV) and 
Boltzmann (7 = 2 meV) results, respectively. Right: Solid, dashed and dash-double dotted lines 
represent non-Markovian {Ti'^ph = 2 meV), Markovlan (7 = 4 meV) and Boltzmann (7 = 4 meV) 
results, respectively. 



B. Electron polarizations 



Figs, m and [5] show the most prominent polarizations between all LL's associated with 
pairs of laser states, obtained from the Markovlan (7 = 2 meV and 7 = 4 meV) and non- 
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Markovian models {h'^ph = 1 meV and h'^ph = 2 meV), as they depend on the magnetic 
field. In all cases, the polarization between the ground state and the upper laser level of 
the next period, shown in Fig. HJ is considerable (~ 10 %), due to the fact that these levels 
actually constitute a doublet state. Although their overlap does not change with magnetic 
field, their polarization does, since the LL electronic structure and all scattering/ dephasing 
processes change as well. The polarizations between the upper laser level or the ground state 
of the previous period and the lower laser level are an order of magnitude smaller (~ 1 %), 
see Fig. |5l This rapidly decreasing trend continues for the polarizations between other pairs 
of states, due to strong electron-LO phonon dephasing. 
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FIG. 4: (Color online). The electron polarization between the ground state of the preceding period 
and the upper laser level (all Landau levels) vs magnetic field. Ns is the total sheet density 
of electrons per period. Solid, dashed, dash-dotted and dash-double dotted lines represent non- 
Markovian (hjph = 1 meV and fijph = 2 meV) and Markovian (7 = 2 meV and 7 = 4 meV) 
results, respectively. 

Generally, nondiagonal contribution curves in each Markovian case are offset to smaller 
values compared to the corresponding non-Markovian ones (7 = 2 meV versus h'jph = 
1 meV; 7 = 4 meV versus h'jph = 2 meV). This is caused by smaller scattering rates 
from populations to polarizations in the Markovian case, since they do not include the 
memory of the interaction. Deviations from this for some magnetic fields are a consequence 
of limited analogy between the two models for the broadening parameters 7 and 7^/^ used 
in the calculation. The coherences in the Markovian case for 7 = 2 meV have larger 
peaks than for 7 = 4 meV, while away from the peak, their values become smaller. The 
effect of a smaller FWHM is that the scattering rates responsible for the formation of 
polarizations have larger peak and lower valley values. The same conclusion applies to 
the non-Markovian results for different damping parameters. To illustrate this, the most 
influential scattering rates from populations to the polarization between the states of the 
doublet in the Markovian case versus magnetic fleld are shown in Fig. [61 The average 
scattering rate from the polarizations between LL's originating from state and rrii^ into 
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FIG. 5: (Color online). The electron polarization between the lower laser level and other QCL 
states (all Landau levels) vs magnetic field. States 1', 2, and 3 represent the ground state of 
the preceding period, the lower laser level, and the upper laser level, respectively. A''^ is the 
total sheet density of electrons per period. Left: Solid and dashed lines represent non-Markovian 
(^7p/i = 1 meV) and Markovian (7 = 2 meV) results, respectively. Right: Solid and dashed lines 
represent non-Markovian {hjph = 2 meV) and Markovian (7 = 4 meV) results, respectively. 



the polarizations between LL's originating from states rrii^ and rrii^ may be defined in the 
Markovian approach according to 
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Under the condition that = mjj, Eq. fl33p gives the average scattering rate from LL's 
associated with state into the polarizations between LL's associated with from states 
rrii^ and m^g. Also, the condition = nii^ and m^g = m^^ gives the average scattering rate 
from LL's stemming from state into LL's stemming from state m^g, which is equivalent 
to the Boltzmann scattering rate. 



C. Optical gain 

The gain spectra for a magnetic field of 4 T in the energy range close to the optical 
transition energies, and one LO phonon energy, are shown in Figs. [7] and [8], respectively. The 
gain was calculated for the non-Markovian {h'-fph = 1 meV and h'jph = 2 meV), Markovian 
and Boltzmann (7 = 2 meV and 7 = 4 meV) dynamics. The modal gain, defined as 
Qmi^^) = g{^)d, and obtained from the Boltzmann theory, has identical major features as 
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FIG. 6: (Color online). Average scattering rates (see text for explanation) from QCL states to 
polarizations between QCL states vs magnetic field, calculated in the Markovian approach for 
7 = 2 meV (solid lines) and 7 = 4 meV (dashed lines). States 1, 2 and 3 represent the ground 
state, the lower laser level and the upper laser level, respectively. States 1' and 3" denote the 
ground state of the preceding period, and the upper laser level of the following period. 



the one predicted from the Markovian approach, for the same value of FWHM. Disregarding 
nondiagonal scattering and dephasing processes in the Boltzmann model, which affect the 
gain linewidth, results only in a quantitative modification of the Markovian prediction. 
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FIG. 7: Optical gain vs energy for a magnetic field of 4 T. The energy range is in the vicinity 
of the optical transition energies. Left: Solid, dashed and dash-double dotted lines represent 
non-Markovian {hjph = 1 meV), Markovian (7 = 2 meV) and Boltzmann (7 = 2 meV) results, 
respectively. Right: Solid, dashed and dash-double dotted lines represent non-Markovian {hjph = 
2 meV), Markovian (7 = 4 meV) and Boltzmann (7 = 4 meV) results, respectively. 

In the case of the non-Markovian dynamics, the gain linewidth is significantly decreased 
for optical transition energies in comparison to the corresponding Markovian and Boltzmann 
estimates (~ 10 times). This result might seem somewhat unexpected because scattering 
and dephasing are increased compared to the Markovian treatment by including the memory 
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of the interaction process. However, in the non-Markovian approach, the broadening caused 
by the scattering terms is energy dependent.— The frequency dependence of the quantity 
^Ki^i^i.^i_^{uj) (see Eq. ( !29|) ). correlated to phonon-assisted matrices, shows that the non- 
Markovian calculation does not yield energy conservation in the scattering processes, but 
includes an additional huj contribution. Then, in the coefficient associated with An|2j)|3.j)(ti;) 
in Eq. (l29l) . if huj E3 — E2, the sum of the terms reciprocal to Ei.^ + hui^o ~ E\2,j) — 
Tioj — ih'jph ~ Ei^ — E\2j) + 23.65 meV—ih'yph and -E|3j) + hu-^o — Ei^ — Tluj — ih'jph ~ 
E\3j) — Ei^ + 23.65 meV— ih'jph appear. For the electronic structure of the QCL considered, 
these terms are not resonant for the majority of = 1,...,2N, and their contribution to 
the coefficient related to Ar;,|2j)|3j)(ti;) is small. Hence, An\2j)\3j){uj) is large, and so is the 
gain. In contrast, the same coefficient, but in the Markovian description, contains the sum 
of the terms proportional to S{Ei^ + Tioji^o ~ and 5(-Eig — ^i^lo ~ -£"14) (see Eq. (I30l) ). 
which are always resonant for 24 = |l,j), or 24 = |3",j), and 25 = |2, j) (or vice versa). 
Therefore, the coefficient associated with An\2j)\3j){uj) is large, and An\2j)\3j){uj) and the 
gain are small. Similar arguments may be used to explain an increased gain peak and 
narrow linewidth in the non-Markovian treatment if hu ~ Ey — E2. However, in a real 
QCL device, it is likely that interaction of electrons with impurities, interface defects and 
other electrons, as well as imperfect periodicity, will broaden LL's making the gain linewidth 
for optical transition energies not as narrow as the non-Markovian model predicts. On the 
other hand, for hu ^ E2 — Ei ^ Tiuji^o in the non-Markovian treatment, the terms inverse to 
Ei^ + hui^o - E\i^j) -hu- ihjph ^ Ei^ - Ell J) - ihjph and E\2,j) + hu^o - Ei^ -hu- ihjph ^ 
E\2,j) — — iTijph contribute to the coefficient related to An\ij)\2,j)iy>) ■ Obviously, these 
terms are resonant for = or ^3 = |2,j), thus An\i j'^\2,j){^) and the gain are small 

and comparable to the Markovian case. The same situation occurs for Tiuj k, E2 — Ey. 

In the Markovian limit, energy renormalizations, describing the polaron corrections to the 
band structure, are ignored. However, the polaron shift is always included in the quantum- 
kinetic treatment. It is more prominent for the energy transitions close to one LO phonon 
energy (~ 1 meV), but it is also present for the optical transition energies (~ 0.4 meV). 

Fig. [9] illustrates the gain profile for a magnetic field of 6 T in the energy range around the 
optical transition energies and one LO phonon energy, calculated from the non-Markovian 
model for h'-fph = 1 meV and h'-fph = 2 meV. Comparison of the non-Markovian gain spectra 
for different values of the damping parameter (see also Figs. [7] and [S]) reveals pronounced 
differences between both the gain linewidth and peak values. Sensitivity of the results to the 
values of phenomenological parameters confirms the need for a self-consistent incorporation 
of higher order correlations in the quantum-kinetic model, which would require a significant 
increase in computational time. Moreover, for a magnetic field of 6 T, apart from the two 
expected peaks associated with the transitions from the lower laser level to the ground state 
or the upper laser level of the subsequent period, an extra peak appears if hjph = 1 meV 
is taken. Careful inspection of Eq. 0291) reveals the presence of more resonances at the 
energies of hu ^ ±{E\yi j'^ — + hui^o = ±2.6 meV+Zic^Lo- At first, it may appear 

that one of these polaron resonances, related to the transition between the ground state 
and the upper laser level of the subsequent period, is manifested via that additional peak 
in the gain/absorption spectra (so-called polaron satellite). However, this is not entirely 
the case here. Despite those resonant terms, the coefficient associated with An|3//j)|ij)(a;) 
is large, since — E^yj) — Tluj is large (see Eq. fl29l) ). At the same time, the coefficients 
associated with An\ij)\2,j){!^) and An\yi ^j)\2,j){l^) are smaller, because E\2,j) — -^|ij> ~ Tiui 
and E\2j) — E^yrj^j — huj are small, and, additionally, they also include the aforementioned 
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FIG. 8: Optical gain vs energy for a magnetic field of 4 T. The energy range is in the vicinity of 
one longitudinal optical phonon energy. Left: Solid, dashed and dash-double dotted lines represent 
non-Markovian {hjph = 1 meV), Markovian (7 = 2 meV) and Boltzmann (7 = 2 meV) results, 
respectively. Right: Solid, dashed and dash-double dotted lines represent non-Markovian {hjph = 
2 meV), Markovian (7 = 4 meV) and Boltzmann (7 = 4 meV) results, respectively. 



polaron resonant terms, as well as the terms reciprocal to Ei.^ + huji^o — —fiuj — ih'^ph ~ 
Ei^ - E\i^j) - ih'jph or E\2,j) + hui^o - - Tiuj - ih^ph ~ E\2j) - Ei^ - ih'jph, resonant for 
^3 = |1, j) or 23 = |2, j). As a consequence, An\ij)\2j){uj) or An\3iij)\2,j){^) constitute much 
larger fraction of the total gain than Ar;,|3"j)|i ,,)(ci;), unlike it is the case in Ref . [isl, where only 
two subbands were considered. Here, the additional peak for B = 6 T and hjph = 1 meV 
is actually related to the transition between the lower laser level and the upper laser level 
of the next period (the calculation showed that An\ij)\2,j){^^) < '-^n\3",j)\2,j){^)), as well as 
the peak at hu = 32.7 meV which is also present for h'-fph = 2 meV, while the peak at 
hu = 36.8 meV is associated to the transition between the lower laser level and the ground 
state. In fact, the previous analysis of the terms which are resonant for the energies in the 
vicinity of one LO-phonon energy suggests that up to 6 peaks could emerge in the absorption 
spectra in that energy range. Due to the their proximity, some of them may be obscured 
by other stronger peaks for different values of the damping parameter and magnetic fields. 
In Fig. [8], for the value of the damping parameter of h'jph = 1 meV, we can notice that, in 
addition to 2 well defined peaks, 4 small ones also appear. 

Fig. [To] shows peak gain values versus magnetic field dependence, for transitions from 
the upper laser level and the ground level of the preceding period to the lower laser level, 
obtained from the non-Markovian {h'jph = 1 meV), Markovian and Boltzmann (7 = 2 meV) 
approach. The general trend of the non-Markovian gain is fairly similar to the Markovian 
or Boltzmann results, although the former has larger values. The oscillations of the gain in 
both types of transitions reproduce reasonably well the oscillations of the related populations 
(see Fig. [3]). It should be noted that, for most magnetic fields, the dominant gain spectral 
component is not generated in the transitions from the upper laser level, but from the ground 
state of the preceding period, even at some magnetic fields for which the upper laser level is 
more populated. The reason for that is a slightly larger dipole matrix element between the 
ground state of the preceding period and the lower laser level. 
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FIG. 9: Optical gain vs energy for a magnetic field of 6 T calculated using the non-Markovian 
approach for ^17^/4 = 1 meV (solid line) and h'^ph = 2 meV (dashed line). Left: The energy range 
is in the vicinity of the optical transition energies. Right: The energy range is in the vicinity of 
one longitudinal optical phonon energy. 




FIG. 10: Maximal gain vs magnetic field for transitions from the upper laser level and the ground 
state of the preceding period to the lower laser level. Left: The non-Markovian {hjph = 1 meV) 
results are represented by solid and dashed lines, respectively. Right: The Markovian (7 = 2 meV) 
results are represented by solid and dashed lines, respectively. The Boltzmann (7 = 2 meV) results 
are represented by dash-dotted and dash-double dotted lines, respectively. 



D. Current 



The current densities as functions of magnetic field, calculated using the non-Markovian 
(^7ph = 1 meV and hjph = 2 meV), Markovian and Boltzmann description (7 = 2 meV 
and 7 = 4 meV), are shown in Fig. [TTJ From Eq. ( l23l) . used in the Markovian and non- 
Markovian approach, it follows that diagonal density matrix elements do not contribute 
to the total current .-^^ Therefore, the electron transport is entirely due to nondiagonal 
density matrix contributions i.e. scattering induced phase coherences between the laser 
states. This quantum-mechanical picture of completely coherent current is in a stark contrast 
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with the semiclassical picture of transport through scattering transitions. However, both 
descriptions give similar results, see Fig. [TTl Here, the nondiagonal density matrix element 
are considerably smaller in comparison to the diagonal ones (see Figs. [3], H] and [5]), thus 
the former may be approximated in terms of the latter,— giving Eq. (1241) used in the 
calculation of the semiclassical current, and resulting in comparable values of the current 
density. 




FIG. 11: Left: Current density vs magnetic field dependence. Left: Solid, dashed and dash- 
dotted lines represent non-Markovian {hjph = 1 meV), Markovian (7 = 2 meV) and Boltzmann 
(7 = 2 meV) results, respectively. Right: Solid, dashed and dash-dotted lines represent non- 
Markovian {hjph = 2 meV), Markovian (7 = 4 meV) and Boltzmann (7 = 4 meV) results, 
respectively. 

In the semiclassical interpretation, the electron transport channel from one doublet di- 
rectly to the subsequent doublet is as important as the channel which additionally involves 
the lower laser level. The scattering rates from the lower laser level to the ground state and 
the upper laser level of the next period do not exhibit pronounced oscillations with mag- 
netic field since those energy transitions are close to one LO phonon energy. Therefore, the 
current versus magnetic field dependence in the semiclassical picture is mainly determined 
by the scattering rates between the doublet states, shown in Fig. [T21 and their populations. 
In the Markovian and non-Markovian description, the current is completely determined by 
the polarizations between the QCL states, see the diagrams shown in Figs. HI and [TTJ In 
both approaches, the current density curves reproduce well the main features of the related 
polarization curves. Also, the discrepancies between the current densities estimated from 
the Markovian and non-Markovian treatments are identical to those between the related 
polarizations. 



IV. CONCLUSION 

We presented a quantum kinetic description of electron dynamics and gain in QCL's 
subjected to a magnetic field, based on the density matrix formalism. As a first step, 
electron-LO phonon interaction was considered as the most relevant scattering/dephasing 
mechanism. Nonequilibrium stationary state populations and polarizations for an example 
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FIG. 12: Average scattering rates (see text for explanation) between QCL states vs magnetic field, 
calculated in the Boltzmann description for 7 = 2 meV. States 1, 2 and 3 represent the ground 
state, the lower laser level and the upper laser level, respectively. States 1' and 3" denote the 
ground state of the preceding period, and the upper laser level of the following period. 

QCL structure were calculated using various kinetic models (non-Markovian, Markovian and 
Boltzmann). Since all of these models contain a phenomenological parameter, we chose the 
sets of those parameters so that steady-state populations are similar for a range of magnetic 
fields, and then we compared other relevant quantities (polarizations, gain, current). In both 
the Markovian and non-Markovian approach, coherent polarizations induced by electron-LO 
phonon interaction were found to be relatively small. This in turn led to similar values of 
the entirely coherent current in the Markovian and non-Markovian picture compared to the 
values of entirely incoherent current in the Boltzmann interpretation. Gain spectra in the 
non-Markovian treatment showed considerably narrow linewidths for optical transitions and 
evidence of polaron formation, in contrast to the Markovian and Boltzmann predictions. 



This appendix presents a detailed derivation of the quantum kinetic equations for electron 
populations and polarizations in QW's under an applied magnetic field. The procedure for 
obtaining appropriate expressions in the Markovian approximation is also given. 

Eqs. f llip and (ITB]) represent a starting point for the derivation of the quantum kinetic 
equations of motion, giving the following expressions 
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Simple algebra then leads to the equations for the dynamics of populations and polarizations 
in the form 
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From the expression for the lateral overlap integrals derived in Appendix B (Eqs. (IB7p 
and (1B9P ). it can be shown that the following expressions 
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where 6 = arg(g^ + iqy) and 



ai 
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hold. Then, the temporal evolution of the density matrix elements read 

«3>«4,«5,q 

-\Hj.,j.,{q.y)\\H,^^,M\Gl^^^^^ ,(A6) 

where the expression for 02 is similar to ai, with the exception that it contains index ji^ 
instead of jj^ . 

Close inspection of Eqs. (lASp and (1A6P shows that they represent independent identical 
subsystems of equations for each wave vector k, therefore the density matrix elements fi^^i^^k 
and the quantities associated to phonon assisted matrices 5K]}^^^^^^^ do not depend on the 
wave vector k under the influence of electron-LO phonon interaction alone. It is obvious 
from Eqs. f|T2l) and f|T3l) that the effect of free electrons system and interaction with light 
on the density matrix elements is identical. Furthermore, since the thermal equilibrium of 
phonon is assumed, phonon populations do not depend on the wave vector q, and so does 
not 6Ki'^i%^ {SKl'^'k% SKi,,,i,i^). Consequently, we have 

r _ _ _ Sfc' fhi2,k' 



eB L^Ly 



d 



«a,«4,«5,q 

{\H,^^,^,^{q,y)\\H,^^,^^{q,y)^^^^ 

dt ih 



5*4 



+ -T [('^o + ^)ni^iA^u,i2 - Oisrii^i^) - noUi^i^iSi^^ir^ - asni^i^J] ,(A7) 
where the substitution asSKi-^^i^i^i^ 5Ki^i^i^i^ is introduced. It follows from Eq. (lASp that 
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if the condition ji^ ^ = ji.^ + jjg is fulfilled, then 0,^ 3 = 0. Furthermore, using the identity 

r2Tr 
=0 



d9 e*f'^i-2e^(^H.2+Ju-i-3-%) = e^5"^-^27r(5,-- +,■- , (A8) 



and summing over the phonon wave vector q, we finally get the expressions given by Eq. f|T6]) 
in Subsec. Ill A[ 

In order to perform the Markovian approximation, it is assumed that the dominant time 
dependence is given by the exponential in Eq. flA2p and therefore the value of electron 
populations can be taken out of integral.— Also, the fast oscillations of polarizations have 
to be taken into account,—"^ resulting in 

^4!^,k~gy = 9*ciH*^j^^{k, k - qy, q^)G*^^^^^^{q,)6{Ei^ + huLo - Ei^) 

[(^0 + '^)fiii5,k{^i4,i2 ~ fui2,k-qy) " nofi4i2,k-qy{Sii,u^ - fhi^,k) (A9) 

Following the identical procedure as in the derivation of the quantum kinetic equations, we 
obtained the equations of motion in the Markov limit given by Eq. flTSl) in Subsec. Ill A[ 

APPENDIX B 

In this appendix we derive the expression for the lateral overlap integral Hj^j-{kf, ki, q^) = 
J u*^ {x — f{kf)) e^'^'^^Ujj {x — f{ki)) dx, with the harmonic oscillator wave function of the 
form 

u,{x) = ^-^HAf3{x - xo)]e-^(— °^ (Bl) 



where P = eB/h, x^ = ki/P"^, Xq/ = kf/P"^, and ki = kf ^ qy (the upper sign holds for 
phonon absorption and the lower one for emission). Inserting t = x — XQf, to = ^^^r^, and 
z = t — to into Eq. (IBip . the lateral overlap integral may be transformed into 



H.Mkf, kf T qy, q.) = ^ , e-^-^o/e-^e^^^ 

j dz H^^{(3z + (3to)Hj^ (^(3z + (3to ± e"^'^'. (B2) 
In order to obtain the final expression for Hjj,j.{kf, ki, g^), we use the following identity 

dx e-"'^''i/™(a + cx)Hnib + cx) = ^ ^' L;^-"^(-2a6), n > m. (B3) 

In the case when ji < jj, by substituting c = m = ji, n = jf, a = [3tQ ± qy/ [3, and 
b = pto, Eq. dBl]) takes the form 

Hj.jAkf, kf T qy, q.) = ... . e-'-^o/e'^e^^^ ^ 



Since the following identities 



2 



2 ' 



hold, the final form of Eq. flB4p reads 
where 



Similarly, if ji > jf, it can be shown that 



^4— / ,.2 \ „2 

xy 



with 

\H,M\ = V-7^] (^) ^: 



1 ^i~^f 



3f3.yH.yn . .,1 I 2/32 I 12/32, 
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